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Abstract 

We calculate the influence of dynamical exchange effects on the response 
properties and the static properties of a two-dimensional many-polaron gas. 
These effects are not manifested in the random-phase approximation which is 
widely used in the analysis of the many-polaron system. Here they are taken 
into account by using a dielectric function derived in the time-dependent 
Hartree-Fock formalism. At weak electron-phonon coupling, we find that 
dynamical exchange effects lead to substantial corrections to the random- 
phase approximation results for the ground state energy, the effective mass, 
and the optical conductivity of the polaron system. Furthermore, we show 
that the reduction of the spectral weight of the optical absorption spectrum 
at frequencies above the longitudinal optical phonon frequency, due to many- 
body effects, is overestimated by the random-phase approximation. 
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I. INTRODUCTION 



Although the static and dynamic properties of a single polaron have been [1,2] (and still 
are [3]) extensively studied, the properties of an interacting polaron system, e.g. their optical 
absorption [4,5], are less well understood. Thus far, the many-polaron system has been 
studied mainly in the Hartree-Fock approximation and in the random-phase approximation 
(RPA) [6-9], even though in Ref. [10] it was shown that in three dimensions (3D) the polaron 
energy and effective mass are influenced by screening effects beyond RPA. In this paper, we 
investigate dynamical exchange effects in both static properties (energy and effective mass) 
and response properties (optical absorption) of the two-dimensional (2D) polaron system. 

However, although our formalism is appropriate for the description of an interacting 
electron (or polaron) gas, it does not allow to investigate the 2D Wigner solid phase formed 
for example by electrons on a helium surface. Also the study of the many small-polaron 
gas and bipolarons in quasi-2D systems such as the copper oxide planes in high-T c cuprate 
superconductors [11-14] lies beyond the scope of our approach in its current form. 

A convenient formalism to take into account many-body effects in the polaron system 
relies on containing these effects in the structure factor [6]. This approach, reminiscent 
of the Feynman-Bijl treatment of superfTuid helium-4, was first developed to calculate the 
ground state energy of a many-polaron system [6], and has later been applied to the optical 
absorption of this system [15]. The 2D dynamic structure factor can be written as 

(1) 




nv q 



with v q = e 2 / (2£oog) the Fourier transform of the bare Coulomb interaction in 2D, n the 
surface density of charge carriers, the dielectric constant at high frequency, and e(q, u) 
the dielectric function of the system of charge carriers. Previous work [6-9,15] relied on the 
RPA dielectric function: 

e RPA (q, u) = l + Q (q, u), (2) 

where Qo(q,u) is the 2D Lindhard polarizability and is known in closed form [16]. A more 
accurate result for the dielectric function is obtained by introducing the frequency dependent 
local field correction G(q, uj) : 

e(q, u = 1 + — r— r. 3 

1 - G{q,u)Q {q,u) 

In the 3D case, the local field correction was derived in the framework of the dynami- 
cal exchange decoupling (DED) method [17]. This method, based on the time- dependent 
Hartree-Fock formalism, was recently extended to the 2D case in Refs. [18,19]. This leads 
to an improved dielectric function, £ded(qs w), that can be used to include dynamical ex- 
change effects in the analysis of the many-polaron system. The DED result for the dielectric 
function is markedly different from the results obtained in the RPA formalism, Hubbard's 
formalism [20] and the formalism of Singwi et al. [21]. In 3D these different approaches can 
be compared to experiment: measurements of the dynamic local field correction by Larson 
and co-workers [22] have shown that the 3D DED approach [17] leads to superior results. 
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In this work, we use the DED dielectric function in combination with the formalism 
describing the many-body effects in a polaron system through the structure factor. In 
section II, we investigate the dynamical exchange effects on the optical absorption of the 
many-polaron system. From the optical conductivity the effective mass of the polaron as a 
function of the density is derived through a sum-rule in section III. Finally, in section IV, 
the ground state energy of the interacting polaron system is calculated beyond RPA. As we 
shall show, including dynamical exchange effects leads to non-negligible corrections to all 
these quantities. 



II. OPTICAL ABSORPTION OF THE INTERACTING 2D POLARON SYSTEM 

The optical absorption of an interacting system of Frohlich polarons in the weak-coupling 
regime was investigated in [7,15]. The optical conductivity in 2D is given by [15] 

9 OO 

Re[(7 2D (^)] = «— J dq q 2 S 2 D(q,uj- uj lo ), (4) 
o 

where e is the electron charge, wlo is the longitudinal optical (LO) phonon frequency and 
a is the dimensionless Frohlich coupling constant determining the strength of the electron- 
phonon interaction. Note that a delta-function peak at uo = 0, omitted in (4), is present 
in the optical absorption. The influence of the Fermi statistics and the screened Coulomb 
interactions between the electrons are taken into account through the structure factor S 2 d 
appearing in expression (4). Thus far, the effect of electron-electron interactions on the 
optical absorption of a many-polaron gas has been investigated only in the framework of 
RPA [7,15,23]. Nevertheless, expression (4) is not restricted to the RPA approximation: any 
form of the dynamic structure factor can be introduced in the integrand. Here, we will use 
the results obtained for the dielectric function in the DED formalism [19] instead. 

Using expression (1) and the DED result for the dielectric function, we can rewrite (4) 

as 

J\f °° / 1 \ 
Re [cr 2D {y)\ = -ot-f—. yT / dk k 3 Im , (5) 

where now k = qjkp and v = hu/(2Ep) are dimensionless wave numbers and frequencies 
based on the Fermi wave vector kp and the Fermi energy Ep, respectively. The prefactor 



^ = n {Tiulo) 2 e 2 / % 47T£oofe| 
2m 2 b ul V 2m b u LO e 2 

contains material constants (the band mass mj,, the LO phonon frequency ujpo and the high- 
frequency permittivity £00) and has dimensions of optical conductivity. Expression (5) for 
the optical absorption allows to easily substitute available results for the dielectric function. 

The DED structure factor can be interpreted as a sum of two contributions, namely a 
contribution from the single-particle excitations (the Landau continuum) and a contribution 
from the plasmon excitations. We can write 



3 



Im 



£ded (k, v) 



-A pl (k)5[v - u p i(k)] + Im 



£cont (k,u) 



(7) 



where z/ L o is the LO phonon frequency in units of 2E F /h and v p \{k) is the plasmon dispersion 
defined by e(k, v p \) = 0. The plasmon excitations in our approach lead to delta-functions, 
with a strength A p i(k), in the frequency dependence of the structure factor. Also the optical 
conductivity can be separated in a contribution due to the single-particle excitations and a 
part coming from the plasmon branch: 

Re[<7 2D (^)] = Re[cr^(v)\ + Re[a c 2 ^(u)], (8) 

with 

X3 r pW m M klA pX (ko) 

Re [o%v{v)\ = a— - , 9 

(z//z/ LO ) 3 v pl (k ) 

Re[a c 2 ^{u)] = -a—. ^ / dk k 3 Im . (10) 

(V/VL0) 3 { \ £ cont (k, V - V LO ) ) 

Here, k is the wave vector at which v — u LO = v p i(k), and v pl = dv p \jdk is the derivative of 
the plasmon frequency with respect to the wave vector. In the 2D electron gas, the plasmon 
branch v p \(k) lies close to the edge of the Landau continuum, k 2 /2 + k. Figure 1 compares 
the plasmon dispersion in the DED approach to the dispersion in the RPA approach for 
several values of r s . For small k values the plasmon frequency in DED is lower than in RPA. 
Furthermore, in DED longer wave length plasmons are present. The oscillator strength 
A p \{k) of the plasmon branch can be obtained straightforwardly using the /-sum rule. If we 
substitute (7) in the /-sum rule, 

1 ^ j 7iTsk di\ 




£T>mi{k,v)) 2^2 
we can find an expression for the oscillator strength of the plasmon branch: 

A » {k) - m {vi - m ) • (i2 » 

where 

°° / 1 \ 

M 1 = - / vlml 71 -^)diy, (13) 

J Q \£cont{k,lS) J 

is the first frequency moment of the Landau continuum contribution only. In Figure 2, we 
compare the oscillator strength of the plasmon branch in RPA and DED for several values 
of r s . At small k the difference between the strength of the plasmons in RPA and DED 
is negligible. As k increases, a difference appears: the maximum oscillator strength for 
the plasmon branch in DED is reached at a smaller wave vector. Compared to RPA, the 
oscillator strength of the DED plasmons shows a long tail as a function of k. 
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In Figure 3 the optical absorption due to the plasmon contribution is shown for several 
relevant values of the density n of a two dimensional polaron gas, using the material pa- 
rameters for GaAs (e^ = 10.8, m b = 0.0657 m e , hu^o = 36.77 meV, a Bohr = 8.78 nm). 
The presence of the shorter wave length plasmons in DED as compared to RPA shifts the 
spectral weight to higher frequencies. Nevertheless, both in RPA and DED the contribution 
of plasmons to the optical absorption of the many-polaron gas is small in comparison with 
the contribution of the single-particle excitations. The total optical absorption is shown in 
Figure 4, again for different surface densities of polarons in GaAs. Although the Landau 
continuum occupies the same region of the {k, z/}-plane in the DED and RPA approaches, 
the spectral weight of the structure factor is distributed differently. As can be seen in Figure 
4, this leads to a measurable difference between the DED and RPA results. Comparing both 
the DED and RPA optical absorption with the single-polaron absorption [24], we find that 
the RPA result overestimates the damping due to many-body effects. 



III. EFFECTIVE MASS OF THE POLARON IN DED 

From Figure 4, it is clear that the spectral weight of the optical absorption in the region 
beyond the single-phonon gap (uj > ujlo) is different for each of the approaches (DED, RPA 
and single-polaron). This may seem surprising at first glance, since the /-sum rule imposes 
that the total spectral weight is a constant: 

7r Ne 2 

Re[a(u)}(kv = ——. (14) 
2m b 
o 

It is important to keep in mind that the optical absorption of the polaron system shows a 
delta- function in the origin (uj — 0), not shown in Figure 4. The total spectral weight is 
distributed between this delta-function peak and the optical absorption spectrum beyond the 
single-phonon gap [25]. The spectral weight of the delta-function is related to the polaron 
effective mass m* in such a way that [25] 

Thus, a change in the spectral weight allotted to the u > ujlo region of the optical con- 
ductivity indicates a change in the effective mass of the polarons. The electron-phonon 
coupling increases the effective mass of a single polaron by Am = m* — m b . The well-known 
single-polaron result for this effective mass increase in 2D is Am/m b = an/8. In Figure 
5, we show the effective mass increase Am/(amj,) = (m* jm b — I) /a as a function of the 
surface density of polarons in GaAs. Both the RPA result (dashed curve, cf. [7-9,26]) and 
the DED result (full curve) approach the one-polaron result in the limit of small density. As 
the surface density of polarons increases, Am decreases, more rapidly in RPA than in DED. 

The fractional change Attided / Am-RPA is shown in the inset of Figure 5, and indicates 
that there is roughly a 15% difference between the effective mass increase predicted by DED 
and that predicted by RPA. This is consistent with the increase in spectral weight of the 
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optical absorption at uo > uj^o in DED as compared to RPA: if less spectral weight is present 
in the region 00 > uj^q, the delta-function in the origin must carry more spectral weight, 
indicating an increase in effective mass. 



IV. GROUND-STATE ENERGY OF THE INTERACTING 2D POLARON 

SYSTEM 

A general expression for the ground state energy per particle of an interacting polaron 
gas at weak coupling was derived in Ref. [6] : 

2m b S 2D (q) 

In this expression E e \ is the energy contribution per particle of the system of charge carriers, 
without electron-phonon coupling, and the second term is the polaron contribution. The 2D 
Frohlich interaction amplitude appearing in this expression is 



l T/ ,2 ihuLo) 2na / h_ 

with A the surface of the 2D system. Whereas the optical conductivity depends on the 
dynamic structure factor, the ground state energy depends on the static structure factor 

°° dw 



S 2 D(q) = [ —S 2 r>(q,u) 

J 7T 





00 



h 



. Im — duo. (18) 

nnv q J \e{q,uj)J 

Expressed in dimensionless variables for the wave number and the frequency, this is 

s ^-~h(^h)) dv - (19) 

A comparison between the static structure factor in RPA and DED is shown in Figure 6. 
Around k = 2, the DED static structure factor shows a transition to the long wave length 
behavior [18,19]. 

The shift in ground state energy per particle due to the electron-phonon coupling AE po \ = 
E(a) - E(a = 0) is given by [6]: 

A^-a^/^ aJggp . (20) 



with 7 = kp/kLo and kio = y^ m b^Lo/^- The result for the polaronic shift in the ground 
state energy per particle, as a function of the Wigner-Seitz radius r s , is shown in Figure 



6 



7 for different approaches. The material parameters for GaAs were used for this figure. 
If no many-body effects are taken into account, the 2D polaronic energy per particle is 
— {it /2)ahu)Lo m the weak coupling regime. If many-body effects are taken into account 
on the level of the Hartree-Fock approximation, we get a monotonous decreasing energy 
as r s increases, while the RPA approximation leads to a minimum in the polaronic energy 
at r s = 1.59, with E/N = —0.70 ahuiLo- We find that including dynamical exchange 
effects lowers |Ai? po i| with respect to RPA, and shifts the energy minimum to slightly lower 
density: r s = 1.61, where E/N = —0.77 aTiujLo- So, using an improved dielectric function 
as compared to RPA leads to a change of ca. 10 % in the ground state energy. 



V. CONCLUSIONS 



Efforts to describe many-body effects in the polaron system have thus far relied strongly 
on the RPA approach [6-9,23,15], even though there have been indications that inclusion of 
a static local field correction modifies the ground state energy and the effective mass of the 
polarons [10]. In this paper, we use a dielectric function that includes a frequency-dependent 
local field correction derived in the time-dependent Hartree-Fock formalism [17-19]. We find 
that taking into account dynamical exchange effects reveals important corrections at weak 
coupling to the optical absorption of the polaron system, the effective mass of the polarons 
and the ground state energy of the polaron system. In GaAs, estimates of the effective 
mass based on RPA are roughly 15% off from the improved results, estimates of the ground 
state energy differ roughly 10%. Furthermore, we show that in the optical absorption, RPA 
considerably overestimates the reduction of spectral weight due to the many-body effects. 
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FIGURES 



FIG. 1. The plasmon dispersion u v \{k) in a 2D electron gas is shown as a function of the wave 
number k (in units of kp), relative to the edge of the Landau continuum k 2 /2 + k, for different 
values of r s . The full curves correspond to the DED result, and the dashed curves to the RPA 
result. In the inset, the location of the plasmon branch and the Landau continuum (hatched area) 
are shown. 

FIG. 2. The oscillator strength A„\ (k) of the plasmon excitation in a 2D electron gas at v p \{k) 
is shown as a function of the wave number k (in units of &f), for different values of r s . The full 
curves correspond to the DED result, the dashed curves to the RPA result. 



FIG. 3. The contribution of the plasmon excitations to the optical absorption of the 
many-polaron gas is shown as a function of the frequency v in units of the LO phonon frequency 
z^LO) f° r r s = 0.64 (corresponding to n = 10 12 cm" 2 ), r s = 0.91 (n = 5x 10 11 cm" 2 ), and r s = 1.17 
(n = 3 x 10 11 cm" 2 ). Full curves represent the DED result, and dashed curves the RPA result. 



FIG. 4. The total optical absorption of a 2D many-polaron system in GaAs is shown as a 
function of the frequency v in units of the LO phonon frequency ^lo, f° r r s = 0.64 (corresponding 
to n = 10 12 cm" 2 ) and r s = 1.17 (n = 3 x 10 11 cm" 2 ). Results are shown for the DED formalism 
(crosses), the RPA formalism (dashed curves) and for the single-polaron case (dotted curves). The 
noise present in the results is due to the numerical treatment. The delta- function in the optical 
absorption at v = is not shown in this figure. 



FIG. 5. The difference between the polaronic effective mass and the band mass Am = m* — nib 
is shown as function of the surface density n of the 2D polaron system in GaAs. The r s values 
corresponding to the densities n are reported on the top axis. The DED result is plotted as a full 
curve, the RPA result as a dashed curve, and the one-polaron result as a dotted line. The inset 
shows the ratio between the effective mass enhancement in DED and RPA, as a function of density. 



FIG. 6. The static structure factor of a 2D electron gas is shown as a function of the wave 
number k (in units of /cp), in the Hartree-Fock approximation (dotted curve), in RPA (dashed 
curves, at r s = 0.5 and 2.0), and in DED (full curves, at r s = 0.5 and 2.0). 



FIG. 7. The polaronic ground state energy in a 2D polaron system in GaAs is plotted as a 
function of r s in the Hartree-Fock approximation (dash-dotted curve), in RPA (dashed curve) and 
in DED (full curve). The dotted line indicates the one-polaron result. 
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